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Red-shifted H 2 emission in NGC3079: More evidence for a 

pc-scale circumnuclear torus? * 
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Abstract. Using the Effelsberg 100-m telescope, sensitive measurements of the H2O megamaser in NGC3079 are 
presented. During 2000 - 2001, 'high velocity' features are seen that are red-shifted up to 225 kms -1 with respect 
to the systemic velocity of the galaxy (Vlsr ~ 1120 kms -1 ). Symmetrically bracketing the systemic velocity, the 
H2O emission covers a velocity range of ~450kms with only one potential narrow gap (~20kms _1 ) near the 
systemic velocity itself. Velocity drifts of individual components are not convincingly detected. It is shown that 
the presence of red-shifted emission and the absence of detectable velocity drifts are not inconsistent with the 
existence of a rotating circumnuclear maser disk at the very center of the galaxy. Significant differences in the 
overall line profile compared to NGC 4258 and a complex morphology of the radio continuum leave, however, 
space for scepticism. 
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1. Introduction 

Recent single-dish and VLBI observations of luminous 
H2O megamasers have been motivated by the discovery of 
an edge-on Keplerian sub-pc scale maser disk enshroud- 
ing a compact supermassive object at the nucleus of the 
LINER galaxy NGC 4258 (e.g. Nakai et al. 1993, 1995; 
Haschick et al. 1994; Greenhill et al. 1995; Miyoshi et al. 
1995; Herrnstein et al. 1999). Circumnuclear disk struc- 
tures traced by H2O are also seen towards NGC 1068, 
NGC 4945, and the Circinus galaxy (e.g. Gallimore et al. 
2001; Greenhill et al. 1997, 2000). H 2 megamasers there- 
fore provide a unique probe to study the kinematics and 
the dynamical structure of the innermost regions of active 
galactic nuclei. 

For the LINER or Seyfert 2 galaxy NGC 3079 (see 
Sawada-Satoh et al. 2000) that contains one of the most lu- 
minous H2O megamasers known to date (e.g. Henkel et al. 
1984; Haschick & Baan 1985; Haschick et al. 1990), argu- 
ments in favor of a circumnuclear disk have been less com- 
pelling. While the known maser components arise in the 
inner few parsecs of the galaxy, they do not align at right 
angles to the radio jet(s), almost all the detected emission 
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* Based on observations with the 100-m telescope of 
the MPIfR (Max-Planck-Institut fiir Radioastronomie) at 
Effelsberg 



is blue-shifted relative to the systemic velocity, and any 
velocity drift of individual H2O features must be small 
(^4.0kms -1 yr -1 ; Nakai et al. 1995; Baan & Haschick 
1996; Trotter et al. 1998; Satoh et al. 1999; Sawada-Satoh 
et al. 2000). 

In this letter we therefore report single-dish observa- 
tions of the 22 GHz water vapor line from NGC 3079, (1) 
to look for the missing red-shifted velocity components 
and (2) to search for velocity drifts in the near systemic 
features that would indicate centripetal acceleration in a 
rotating circumnuclear disk. 



2. Observations 

Observations of the JK a K c = 616 — 5 2 3 H2O maser line (rest 
frequency: 22.23508 GHz) were made with the MPIfR 100- 
m radio telescope at Effelsberg between February 1994 
and December 2001. Until 1998, a single channel K-band 
maser receiver was employed in a position switching mode 
with a system temperature of T sys ~ 225 K on a main 
beam brightness temperature (T m b) scale. An autocorre- 
lator provided bandwidths of 50, 25 or 12.5 MHz with 1024 
spectral channels, yielding channel spacings of 0.66, 0.33 
or 0.16kms~ . Since 2000, we used a dual channel K-band 
HEMT receiver in a dual beam switching mode with a 
beam throw of 2' and a switching frequency of 1 Hz. After 
averaging the two orthogonally polarized signals, T sys ~ 
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Table 1. Effclsberg 100-m observations log a ) 
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a) Several tunings and bandwidths were employed during most 
epochs. As a consequence, noise levels and channel spacings are 
not uniform across the observed velocity ranges, la noise levels 
vary from 200 mjy (channel spacing 1.32 kins -1 ) in September 

1996 to 7mJy (channel spacing 1.05kms _1 ) in December 
2001. Velocity resolution ~ 1.25 x channel spacing. For direct 
access to 64 flux calibrated spectra, contact chenkel@mpifr- 
bonn.mpg.de. 

180 K. Four autocorrelator backends were employed for 
each receiver channel. Each spectrum had a bandwidth of 
40 MHz and 512 channels, yielding a channel spacing of 
1.05 km s -1 . 

Amplitude calibration was based on measurements of 
the 22 GHz continuum flux of 3C 286 (see Baars et al. 1977; 
Ott et al. 1994). Pointing measurements toward nearby 
sources (in most cases DA 251) were made once per hour. 
The resultant pointing accuracy was < 8" , that should be 
compared with the full width to half power beam size of 
40". Calibration uncertainties are estimated to be ±15%. 

3. Results 

H 2 spectra were taken during 16 epochs (Table 1; no 
changes in the maser profiles were seen within any of these 
observing periods). Several H2O features could be traced 
over the entire monitoring period. The most prominent 
component, detected at ^956 km s -1 , showed a peak flux 
density of 2.5-4.5 Jy. While the peak velocity is observed 
to drift from about 955.0 to 956.8 km s" 1 (estimated er- 
ror of individual measurements: 0.3 km s -1 ), this drift is 
not systematic and appears to be caused by a multitude 
of individual, variable subcomponents. During times with 
lower peak velocity, the profile tends to show a 'shoul- 
der' at the high velocity wing of the line (e.g., on Nov. 
5, 1997 and Feb. 1, 1998) and vice versa (Dec. 21, 2000, 
and Apr. 22, 2001). Notable were two flares of more red- 
shifted components that exceeded 2Jy: The 979 km s -1 
component with 0.7 Jy during Feb. 1994 to Apr. 1997 (see 
Nakai et al. 1995 for earlier spectra) reached 0.9 Jy in Nov. 

1997, 1.5 Jy in Feb. 1998, and 2.2 Jy in May 1998 (Fig. la); 
it then decreased to 1.8 Jy in June 1998, 1.25 Jy in Aug. 

1998, and 0.35Jy since Oct. 2000. The 1017kms _1 ve- 
locity component, observed between Sep. 1996 and Nov. 

1997 at 0.1-0.5 Jy, reached a peak flux density of 2.5 Jy 
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Fig. 1. The (a) 979 and (b) 1017kms" 1 features flar- 
ing beyond the 2Jy level (channel spacings: 0.66 and 
1.05kms _1 ). Velocity scale throughout the paper: Local 
Standard of Rest (LSR), following the radio astronomical 
definition of velocity (see e.g. Trotter et al. 1998). Vlsr = 
Vhel + 3.23kms _1 . 

in Mar. 2000 (Fig. lb); flux densities decreased to 1.0 Jy 
in Oct. 2000 and to 0.1 Jy in May 2001. Emission near 
the systemic velocity of the galaxy (1080 km s -1 < Vlsr 
< 1160 km s -1 ) remained faint and not a single narrow 
velocity component could be traced during several con- 
secutive observing epochs. 

Fig. 2 shows a spectrum of near systemic 
(Vsys-mOkms -1 ) and red-shifted {V>V sys ) fea- 
tures that were detected since March 2000. More than 10 
distinct narrow components appear at velocities in excess 
of 1100 km s -1 (Table 2). We conclude that the H 2 
maser emission extends over a radial velocity range of 
~450kms _1 . The integrated luminosity of the previously 
known blue-shifted components remains approximately 
constant. This may also hold for the systemic and 
red-shifted features shown in Fig. 2. A detection of 
these prior to 2000 would have been difficult in view of 
technical improvements that occurred at Effelsberg in 
1999. Sensitivity also limits the studies by Nakai et al. 
(1995; we do not confirm the weak 764 and 791 km s -1 
components) and Baan and Haschick (1996). However, 
the 1123 and 1190 km s -1 components observed by 
Trotter et al. (1998) in Jan. 1995, the flaring 1192 km s -1 
feature (~0.4Jy) detected by Nakai et al. (1995) in 
Apr.-May 1995, and a 1201 kms -1 feature detected by 
us in Sept. 1995 (Fig. 3) indicate the presence of isolated 
red-shifted maser features as early as half a dozen years 
ago. 

4. Discussion 

In NGC4258, Nakai et al. (1993) detected red- and 
blue-shifted H 2 'satellites' symmetrically displaced by 
approximately ±900 km s -1 from the systemic velocity. 
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Table 2. Velocities of distinct near systemic or red-shifted 
features observed on December 5, 2001 (see Fig. 2) 



^(kms- 1 ) 



1062. 
1066. 
1076. 
1123. 
1137. 
1182. 
1192. 



7±0.4 
8±1.4 
5±0.2 
0±0.9 
7±1.4 
6±0.4 
9±0.5 



V^kms- 1 ) 

1209.3± 0.3 
1229.8± 0.4 
1251. 1± 0.1 
1262. 1± 0.2 
1291.5± 0.9 
1314.4± 1.0 
1347.6± 0.4 
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have seen, for the first time, all the stronger 22 GHz H2O 
maser components in NGC 3079. A nuclear systemic ve- 
locity of 1230 kms" 1 as suggested by Satoh et al. (1999) 
is not supported by our H2O data. 

Nevertheless, the maser emission is far from perfectly 
symmetric: (1) All bright components (>0.5Jy) are blue- 
shifted with respect to the systemic velocity which is in- 
consistent with the situation in NGC 4258 and NGC 1068 
as well as with the spiral shock model proposed by Maoz 
& McKee (1998). (2) The systemic velocity obtained from 
model fits (Irwin & Seaquist 1991) lies at the edge of 
that narrow velocity interval (Ay~20kms _1 , see Fig. 2) 
that appears to be devoid of H2O emission (the only such 
interval over the entire H2O velocity range observed by 
us) . Interpreted in terms of the paradigm established for 
NGC 4258, this may hint at a lack of nuclear 22 GHz ra- 
dio continuum emission at the very center of the putative 
masering torus. (3) In a well ordered edge-on circumnu- 
clear disk, small line-of-sight velocity gradients that stim- 
ulate maser emission should only be present near the tan- 
gential points and near the front and back side of the disk. 
A continuous coverage of the velocity range as observed in 
NGC 3079 is not expected. (4) None of the four systemic 
velocities quoted above provides a particularly high num- 
ber of red- and blue-shifted components that match each 
other with respect to velocity (|K e d — Kiysl ~ I ^biue — Vsys | ) ■ 
Note that the following discussion, does not sensitively de- 
pend on the exact choice of V^ ys as long as it is in the range 
1100-1150 km s" 1 . 

Most blue-shifted 22 GHz H 2 velocity components 
of NGC 3079 are already known to show no significant 
velocity drift (e.g. Fig. 7 of Nakai et al. 1995; Fig. 6 of 
Baan & Haschick 1996) . The slow drift of -0.4 km s _1 yr" 1 
suggested by Baan & Haschick (1996) for the 945, 951, 
and 1015 kms -1 components (V refers to epoch 1994.0) 
is not expected if a circumnuclear torus is present (e.g. 
Miyoshi et al. 1995). Our data for the 945 km s -1 com- 
ponent (945.7 and 946.2 km s _1 in Feb. 1994 and 1998, 
and 946.6 and 946.5 km s -1 in May and June 1998; er- 
rors derived from Gaussian fits are ~0.3 kms) are not 
contradicting Baan & Haschick (1996), although a group 
of features with varying amplitudes could also simulate 
the drift. The 1017kms _1 component (Fig. lb) fluctu- 
systemic velocity of — 1120kms _i and suggests that we ates between 1014 and 1018kms _1 and does not sup- 



Fig. 2. Near systemic and red-shifted H2O maser features 
toward NGC 3079. Channel spacing: 2.1 kms -1 (main 
spectrum), 4.2 km s^ 1 (inserts (b) and (c)), and 8.4 km s -1 
(insert (a)). The inserted spectra (insert (c) includes a few 
bad channels at the edge of the band) show an initially 
dominant fading 1222 kms -1 component, while the flux of 
the 1192 kms -1 component increases. The brightest com- 
ponents have an isotropic luminosity of —0.5 L©. Since the 
zero level is determined by a linear baseline derived from 
velocities V < 920 and > 1360 kms - , the absence of H 2 
emission near 1100 kms -1 is not entirely certain. 



This suggested the presence of a circumnuclear disk that 
was later confirmed (see Sect. 1). Is the discovery of a 
multitude of near systemic and red-shifted features in 
NGC 3079 also hinting at a circumnuclear disk? Emission 
is seen at least between 928 and 1352 kms -1 (Fig. 3). 
Assuming the presence of a symmetric edge-on circum- 
nuclear disk, the systemic velocity of the nuclear region 
then is ~1140kms~ 1 . Including the 896kms~ 1 feature 
apparent in the early spectra of Henkel et al. (1984) and 
Haschick & Baan (1985) and tentatively also seen in Dec. 
2000 and Dec. 2001, V sys ~ 1124 kms" 1 . Is this consistent 
with the systemic velocity of the galaxy? 

From the nuclear molecular gas traced by the CO J 
= 1-0 line, Irwin & Sofue (1992) deduced a systemic ve- 
locity of V sys — 1050±10kms _1 . From an Hi line profile, 
Irwin & Seaquist (1991) obtained a midpoint at the 20% 
peak level of 1123±10kms _1 and 1126 kms -1 from a fit to 
the moment map. Three dimensional modeling of their H 1 
data cube yields V sys = 1116kms _1 . The rough agreement 
of CO, Hi, and H 2 velocity centroids implies a nuclear 
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Fig. 3. Individual velocity components of the 22 GHz H2O 
spectrum obtained with an accuracy Sjlkms -1 for the 
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port a regular drift, although it is likely related to the 
1015 km s -1 component of Baan and Haschick (1995). The 
951 km s -1 component was not seen by us. Satoh ct al. 
(1999) and Sawada-Satoh et al. (2000) report a drift of 
3.7±0.6kms _1 yr -1 for a maser feature near 1190 km s -1 . 
We find a group of features (total width: I4 ot ~5-6 km s _1 ) 
centered at 1192.9kms~ 1 , both on Mar 3. and Dec. 5, 
2001. Thus the drift appears to be well below our esti- 
mated detection level of 1.5kms _1 yr _1 for a time inter- 
val of 9 months. In view of the long lifetime inferred for 
some of the blue-shifted maser lines (Sect. 3 and Fig. 3), 
the 1192kms _1 feature reported by Nakai et al. (1995) is 
likely part of this component, reducing the velocity drift 
to well below 1.0 km s -1 yr -1 . 

The overall lineshape of the H 2 profile is quite dis- 
tinct from that of NGC 4258 and the morphology of the 
radio continuum (T. Krichbaum, priv. comm.) is more 
complex. In spite of the presence of a roughly linear ridge 
of blue-shifted H2O masers, the existence of an associated 
nuclear disk is thus not certain (see e.g. Trotter et al. 1998 



for alternative scenarios). Sawada-Satoh et al. (2000) pro- 
posed a thick maser disk with continuum source 'B' at the 
nucleus and the 1190 kms -1 feature on its near side (their 
Fig. 6). Since we see no velocity drift near 1190 kms -1 
and since this is not the systemic velocity, however, there 
is no need to put this component at the near side of the 
putative disk. We therefore tend to favor the disk sce- 
nario outlined by Trotter et al. (1998) with the nucleus of 
the galaxy being located between continuum components 
A and B (their Fig. 7). In this latter scenario, expected 
velocity drifts would be consistent with our upper lim- 
its. The nucleus would contain a few 10 6 M Q within R 
~ 10 mas (~0.7pc; assumed rotation velocity V^ot^lOO 
225 km s -1 ), the centripetal acceleration and proper mo- 
tion of the near-systemic components would be difficult to 
detect (<0.1 kms -1 yr -1 and ^l^asyr -1 ), and most of 
the red-shifted masers would be located 10-15 mas south 
of the blue-shifted ones (some masers may be associated 
with the jet (e.g. Trotter et al. 1998)). Sensitive VLBI 
data are needed to examine this picture and to obtain 
definite information on the relative location of blue- and 
red-shifted maser features and continuum sources in the 
very center of NGC 3079. 
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